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I. INTRODUCTION 

In condensed matter physics, strongly correlated elec- 
trons underly a great variety of ordered states, both com- 
mon and exotic (e.g. ferromagnets, superconductors). 
One of the lesser-studied orders is spontaneous currents 
(known sometimes as "orbital antiferromagnetism"). In 
this paper, we seek a minimal (spinless) toy model that 
manifestly exhibits such currents, precisely because any 
systematic study of fermion orderings with a quadratic 
order parameter reveals that the possible ordered states 
include not only the familiar cases of charge or spin den- 
sity waves or superconductivity, but also spontaneous 
orbital currentsii^ii^ii^i^. Yet such states have not been 
definitively observed in any material, nor numerically in 
the Hubbard modeJS-, and only very recently for any real- 
istic microscopic Hamiltonian^'-. Thus, we ask : Which 
aspects of the interactions and/or degrees of freedom dis- 
pose a system generically towards ordered states with 
spontaneous currents ? 

Such states were considered especially in the context 
of high- Tc cuprates. Early in their history, "flux phases" 
with current order were invente d^'^^d^ however the ac- 
tual phase was expected to be disordered. More recently, 
two different kinds of spontaneous-current order were ad- 
vanced to explain the mysterious pseudogap state of high- 
Tc cuprates^i^. Ref. [sj proposed the "d-density wave", 
which breaks translational symmetry (currents circulate 
in opposite senses around even and odd plaquettes); vari- 
ants were considered more recentljsiS, e.g. modulated 
versionsi^iii. 

In contrast, Varma's phase o'^'^^ require the so-called 
"three band" model in which oxygen orbitals of the Cu02 
layer are explicit independent degrees of freedom; the lat- 
ter state breaks 4-fold rotational and time-reversal sym- 
metry, but not translational symmetry. Experiments on 
photoemission^*' (in BSCCO) and neutron diffraction^^ 
(in YBCO) indicated time-reversal symmetry breaking, 
in the pattern of Ref. [l^ Finally, Khomski and collabo- 
rators showed currents are implied by non-coplanar spin 
order in (spinfuU!) Mott insulatorsi^. 

These proposals motivate a basic question: under what 



circumstances, in principle, can a quantum state be re- 
alized with spontaneous currents? Where, in a model's 
parameter space, is such a state favored? Ever since the 
Hubbard model, toy lattice models having a minimal pa- 
rameter space (and possibly amenable to solution) have 
been key tools to sort out basic questions such as these. 
For the more familiar orders, "strong-coupling" models 
are well-known in which some "zero-order" state trivially 
has the order in question, and the order is stable against 
small perturbations. Thus, in the phase diagram, one 
is assured of a corner where the ordered phase occurs 
and extends an indetermined distance towards the regime 
where perturbations are large (which is usually the phys- 
ical regime). But in the case of currents order, no general 
intuitive picture has emerged. 

This paper addresses this question using a toy-model 
built from square plaquettes; focusing mainly on the sim- 
plest case of spinless fermions, we explore the possibilities 
for realizing spontaneous currents. The main prior study 
of orbital currents in spinless models is Nersesyan's lad- 
der modeli^i22,, in which a map to spinfuU chains was 
introduced that we adopt in Sec. IIVBI Quite recently, 
spinless models were motivated by the possible realiza- 
tion in cold dilute atoms^i. 

The choice of square plaquettes is a choice motivated 
both by convenience of calculation and real material ge- 
ometries. As we will see, a square plaquette has sponta- 
neous currents as one of its natural degrees of freedom 
which is what we desire to investigate: possibility of spon- 
taneous currents in the zero-order ground state. 

This paper is organized as follows: In Section |TT1 we 
define our toy-model Hamiltonian and set up the vari- 
ous lattice geometries - tube, ladder, and square lattice 
- we shall deal with; we go on to describe the properties 
of one square plaquette as it forms the basic unit of all 
the lattice geometries considered, in particular reducing 
its degrees of freedom to a pseudospin via the method of 
canonical transformations (which is briefly summarised 
in the appendix, as it is the basis of all our subsequent cal- 
culations.) The core section is Sec. IIII[ where we imple- 
ment the pseudospin projection (illustrating it in detail 
for the case of a "tube" lattice) and obtaining a pseudo- 
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spin effective Hamiltonian, showing its final form for the 
respective lattices; we also explore the relation between 
the fermion Hamiltonian and the pseudospin Hamilto- 
nian, focusing on possibility of spontaneous currents in 
the ground state. In Sec. llVi we connect out work to 
spinfuU models in two ways: simply incorporating spin 
fSec. lIV A| or mapping a pair of site indices to spin labels 
(Sec. IIVB[) . At last in Sec. |V] we ask if we have learned 
how to construct a uniform lattice model with currents. 
We conclude fSec. lVIl) by discussing why it is hard to ob- 
tain spontaneous current order, and what light this may 
shed on realistic motivated models of such order. 



II. MICROSCOPIC MODEL AND PSEUDOSPIN 
MAPPING 

Our basic model Hamiltonian is H = H^'°p + +H' 
with 

T^hop = -t ^ [c^(r)c(r') + /i.c.] (2.1a) 

n.n.D 

= +y^n(r)n(r') (2.1b) 

□ 

H' = -t' [c^(r)c(r') + /i.c.]. (2.1c) 

n.n.D-n 

Each site r has an orbital with room for one spinless 
fermion. A disjoint set of plaquettes ("strong plaque- 
ttes") are singled out. Within each strong plaquette 
(tagged by "□" in notations) there is a hopping —t on 
every bond; there is also a repulsion V between any two 
fermions (whether first or second neighbors; in a spinless 
model, of course, there can be no onsite term). 

Finally, every bond between bold plaquettes has a hop- 
ping —t' , which is assumed to be a small perturbation. 
We will consider three kinds of lattice geometries, as 
shown in Fig. [1] (a) a ladder, in which every other pla- 
quette is strong; (b) a "tube", which is one-dimensional 
like the ladder, but the strong plaquettes are oriented 
transverse; and (c) a square lattice, in which one of four 
plaquettes is strong. The ladder is simplest, but also has 
the least symmetry. 

a) ^ c) 




b) 



FIG. 1: Geometry of a (a) Ladder, (b) Square Lattice and 
(c) Tube. The fermions reside on the vertices of the lattices 
shown. 

We were guided by three considerations while concoct- 
ing the model of (|2.ip . Firstly, we desired the sponta- 



neous currents to be explicitly related to the degrees of 
freedom describing our (degenerate) ground state sub- 
space. (See subsection III Bl below.) Secondly, the (zero 
temperature) behavior should be obvious in a strong cou- 
pling limit. A standard trick^^"^'^. to achieve both ends 
is to artificially weaken some bonds thereby introduc- 
ing a small parameter (t' in our case). In the t' = 
limit, the system decomposes into small disjoint clusters, 
each with a degenerate ground state whose operators are 
represented by pseudospins;. As the small parameter is 
perturbatively turned on, it generates an effective Hamil- 
tonian between the pseudospins; from the symmetry of 
the effective Hamiltonian, one can often read off the sym- 
metry of its ground state. 

Finally, to make our model more physical, we limit the 
terms to fermion hoppings and interactions and no other 
four-fermion terms. Also, as we hope that our model(s) 
might later be adiabatically connected to a uniform one 
(see Sec El , if a certain term is included (say) within 
strong plaquettes, we will be open to including an inter- 
plaquette term of the same form (with arbitrary - small - 
coefficients). But we never assume any particular condi- 
tion on the ratios between the intra- and inter-plaquette 
terms, except that all of the latter are small for pertur- 
bation purposes. 

A. Eigenstates of Disconnected Plaquettes 

Let 7^0 include the 7i''°P and Ti,^ terms, representing 
a set of disconnected squares. We will work at half fill- 
ing, i.e. two fermions per square on average, but our 
Hilbert space includes all ways of distributing these over 
the plaquettes. 

Consider an isolated strong plaquette, with sites x = 

0. 1.2.3 forming a ring. Note TiX a is the same for all 
states accessible by hopping, so if there are fermions 
on the plaquette, TiXu = \n^{n^ — 1)T^ drops out like a 
c-number: as in a noninteracting model^i, multi-fermion 
states are built from the one-particle eigenstates on the 
ring, defined by creation operators 

4^i^e5™ct(x) (2.2) 

X 

where m — 0, ±1, 2 is the angular momentum around the 
ring. The single- fermion eigenenergies are 

E„, = -2tcos{\TT'm), (2.3) 

1. e., Eq ^ -2t, E±i ^ 0,^2 = +2i. Table U lists the 
multi-fermion ground states for each occupation sector 
of a single plaquette. Our interest will be the 2-fermion 
sector since it has degenerate ground states |2-|-) and |2— ) 
with spontaneous current in the + and — senses, respec- 
tively. 

To have any possibility of a symmetry broken state, (at 
least some of) the plaquettes must be in the degenerate 
half-filled ground states. What is the ground state of 
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TABLE L States with n fermions on a plaquette. 



rP 


label occupation 


energy 










1 


|i> |6> 


-2t 


2 


|2+) |o,+i> 


-2t + V 




|2-> |6,-i> 


-2t + V 


3 


|3> |o,+i,-i) 


-2t + 3V 


4 


|4> |o,+i,-i,2) 


6V 



an extended system of N sites forming N/A disconnected 
plaquettes with N/2 fermions (i.e. half filling)? The case 
y = is more degenerate than we wished, since any 
combination of states with nP — 1,2,3 has total energy 
—2t{N/A). However, taking V > favors the subspace 
in which ri^ = 2 on every plaquette. In that case, the 
only freedom is the senses of the currents in each of the 
N/4: plaquettes, giving a degeneracy 2^/^. 



B. Pseudospin Mapping. 

These states can be labeled as an array of spin-i pseu- 

dospins Pa with = ±1/2 when plaquette a is in state 
|2±). We aim, via second-order perturbation in t' , to 
compute the effective Hamiltonian Ti.^ defined within the 
ground state manifold (and thus taking the form of a spin 
Hamiltonian in {Pa}.) 

The spin-i pseudospin Hilbert space can be defined as 
follows : 

|±), ^ |±) ^ |2±) 

|±). ^ -i=(|2±)±|2T)) 

\±)y ^ i=(|2±)±z|2T)) 

(2.4) 

Different orders of the fermions - spontaneous cur- 
rents, and (site- or bond-centered) charge density waves - 
correspond to expectations of three characteristic opera- 
tors; when projected to the pseudospin subspace (|2.4p . 
these reduce to the three pseudospin operators (here 

1, j = 0, ...,3 label sites counterclockwise around a pla- 
quette, as in Figs. [3] and [4|): 

1. Pseudocurrent operator: 

pz 

4 = = «(4cj - c]c,) — (2.5a) 

2. Charge Density operator: 

n., ^ cU - (-1)'^ + \ (2.5b) 



3. Bond Density operator: 



4- ^ (clc, + c]cO ^ + - (2.5c) 

Here "— >" means the operators have the same matrix 
elements when acting in the pseudospin Hilbert space. 
Any operators in the pseudospin subspace of a plaquette 
can be expressed in terms of P = (P^,P^,P^). Fig. 
[2] depicts states in which the respective operators have 
expectations. 



l2+>z 
1/2 1/2 

1/2 ^ 1/2 


l2+> 

X 

3/4 1/4 
1/4 3/4 


l2+>y 

1/2 1/2 
1/2 1/2 


1/2 \^ 111 

1/2 ^ 1/2 
l2->z 


1/4 3/4 


1/2 1/2 






3/4 1/4 
l2->x 


1/2 1/2 
l2-> 

y 



FIG. 2: Properties of single plaquette eigenstates. Arrows 
represent orbital currents. Lines represent the bond density; 
non-dashed line represents unit bond density, dashed line rep- 
resents half and absence of line represents zero bond density. 
The numbers at the corners of the plaquette represent the 
charge density at the respective sites. Note that they add up 
to two corresponding to half-filling. 

So, for purposes of nomenclature, we call the pseu- 
dospin states in the z-direction as current carrying states 
(which is expected since -|-1 and —1 carry momentum) or 
CCS, the pseudospin states in the x-direction as charge 
density waves or CDW, and the pseudospin states in the 
y-direction as bond density waves or BDW. Spontaneous 
currents. Orbital currents or just currents will be used in- 
terchangeably to refer to CCS as they have been used in 
the literature before. (The "bond order" , making differ- 
ent directions inequivalent without a translational mod- 
ulation of the charge density, would be an example of 
"electron nematic"— if all similarly oriented bonds had 
the same order parameter.) 

Incidentally, we use the term pseudocurrent operator 
because this is not the true current operator. The latter 
would be the time derivative of the charge density op- 
erator, and is evaluated as a commutator of the charge 
density operator with the full Hamiltonian. Most often, 
the true current is proportional to (or at least has over- 
lap with) the pseudocurrent operator; then, any state 
with pseudocurrent order will also have true current or- 
der. (The pseudocurrent would be the real current if the 
Hamiltonian contained only nearest-neighbor hopping.) 
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III. EFFECTIVE PSEUDOSPIN HAMILTONIAN 



Effective Hamiltonian for the Tube 



In this section, we go on to calculate an effective Hamil- 
tonian by second-order perturbation theory, formulated 
via canonical transformations (reviewed briefly in Ap- 
pendix [X]). We shall consider several variations on the 
model Hamiltonian (|2.ip : the result is always a special 
case of the general form 



The perturbation (t' hopping) changes the filling on 
two plaquettes, hence no first order process stays in the 
reduced Hilbert space (of n^a = 2 on all plaquettes). 
To do that in a second order process, a fermion hops 
from one plaquette ("A") to a neighboring one ("B"), 
and then a fermion hops back from the second to the 
first plaquette. 



(a/3> 



Here Pa is the pseudospin; (a runs over all strong plaque- 
ttes and (a/3) are nearest neighbors (in the ladder, tube, 
or square lattice arrangements). In view of the mapping 
of operators (|2.5p . the original system has spontaneous- 
current order if and only if (j3.ip has pseudospin order 
in the z direction. Consequently, our central concern is 
whether and how the terms in ()3.1|) break pseudospin 
rotation symmetry. 

The above form of the effective Hamiltonian is gov- 
erned by the two different kinds of inter-plaquette hop- 
ping processes that can occur at second-order perturba- 
tion theory. As we will see in the next subsection, at sec- 
ond order, only two adjoining plaquettes can take part in 
the hopping processes. On the one hand, a "degenerate" 
hop takes a fermion from a plaquette to a single-fermion 
orbital on the adjoining plaquette, degenerate with the 
orbital it hopped out of; such hops are responsible for all 
the pseudospin exchange interactions. 

On the other hand, an "excited" hop takes the fermion 
from a partly-occupied degenerate orbital on one plaque- 
tte to a higher orbital on the adjoining plaquette, which 
would not be occupied in any single plaquette state; this 
kind also includes the case (related by particle-hole sym- 
metry) of an electron hopping out of a deeper orbital into 
a degenerate orbital on the new plaquette, or even into a 
higher one (and of course the reverse hops). An excited- 
state hop is not conditioned on which state the plaquette 
(with the excited state) was in: within the pseudospin 
manifold, that only differs in the "degenerate" partly- 
occupied levels. Morever, the only way to return to the 
pseudospin manifold is to undo the same hop thus mak- 
ing exchange of pseudospins impossible. Consequently, 
excited-state hops can (at most) generate only single- 
pseudospin terms in the effective Hamiltonian. 

We shall first consider the tube model case (Sec. IIII Ap , 
since it has the greatest symmetry (the combination of 
two adjoining plaquettes has a 4- fold rotation). The 
other two cases (Sec. IIII Bp are variations on the tube 
case, in that either additional terms appear (due to re- 
duced symmetry) or are accidentally canceled. 




FIG. 3: Tfie two-plaquette unit involved in a second order 
hopping process. Site labels show our convention for the (bi- 
partite) ladder or square lattices. 



For a pair of plaquettes on a tube (Fig. [3]) , the pertur- 
bation takes the form 

i 

= -^'(E4™CB.™ + /».c.) (3.2) 



Notice that the hopping conserves the angular momen- 
tum around the plaquette. For this reason, the only 
excited states that can participate in the second or- 
der hopping processes are |1;3) = |0; 0, -1-1, — 1) and 
|3;1) = |0,-|-1,— 1;0) and the corresponding non-zero ma- 
trix elements are 



(1;3|W'|2- 
(3;1|W'|2- 
(1;3|H'|2- 
(3;1|7^'|2- 



2-> 

2+) 
2+) 



t' 
t' 



-t' 
-t' 



The rest of the matrix elements are zero. Thus, using 
Eq. lAll we get the following second order effective two- 
plaquette Hamiltonian 

Htube =-— (|2+; 2-)(2+; 2 - I + |2-; 2+) (2-; 2 4 

-H— (|2+;2-)(2-;2+| + |2-;2+>(2+;2-|j 

(3.3) 

Conversion to spin Hamiltonian: In accord with our 
pseudospin mapping, we abbreviate |2±) by |±) to label 
the pseudospin states. Now, the transcription to spin 
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notation (for pseudospin(P)) is: 
I2+X2+I - (i 
I2-X2-I - (i 
|2+)(2-| 
|2-)(2+| 



P^); 



P-. 



(3.4) 

Inserting Eq. 13.41 into Eq. 13. 3( we get for the infinite 
tube 

4i'2 n 



Hi 



tube 



V 



1 



+ Pa ■ Pa+1 



(3.5) 



Thus, the effective pseudospin Hamiltonian for the tube 
is a one dimensional spin- 1/2 Heisenberg Antifcrromag- 
net which, as is well known, does not exhibit long range 
order but only power-law correlations. 

This calculation is not only reminiscent of, but com- 
pletely analogous to, the derivation of the effective 
Heisenberg antiferromagnetic exchange interaction in a 
half-filled Hubbard model; the role of spin is taken by 
our angular momentum, since it is conserved by the hop- 
ping along the tube. Hence only our | ± 1) single-particle 
states (analogous to spin up and spin down electrons) 
take part in the second-order process, thus giving rise to 
effective pseudospin exchange of exactly the same (rota- 
tionally symmetric) form as spin exchange in the Hub- 
bard model. This is exactly the content of the discussion 
on "degenerate" hops in Sec. IHII 

In other cases of our model (ladder or square lattice), 
the perturbation need not conserve angular momentum, 
so excited states like |0; 0, -|-1, 2) or | ± 1; 0, — 1) may 
then mediate second-order processes via the "excited- 
state" hops defined at the beginning of this section Thus, 
they will give rise to single-site pseudospin terms only. 



B. Ladder and Square Lattice 




FIG. 4: The two-plaquette unit involved in a second order 
hopping process for the ladder. Site labels are shown, as used 
in e.g. 112} . 

In this sub-section, we list down the results of similar 
calculations for Ladder and Square Lattice cases. For the 
Ladder, the effective pseudospin Hamiltonian is 

+'2 



dder 



y— 

^ V 

a. 



(3.6) 



where x = 2t/V . To make the symmetry of the above 
expression clear, we make a simple transformation as fol- 
lows. 

Staggered Pseudospins: Let us define a new set of 
staggered spin operators Tq , by switching the definitions 
of "up" and "down" pseudospin on every other site "B" 
by a 180° rotation around y-axis. Then, 



py 
P± 



Bi 



^B^ 
^ B 



(3.7) 



while pseudospin operators on sites "A" stay unchanged. 
This transformation converts Eq. (|3.6p to 



H„ = 



E 



— I const[x) 



- Ta ■ Pa+i + , 



2V{1 



(3.8) 

The need for staggering the pseudospin arose out of a 
technicality. Let's focus on the effect of one pair of hops 
connecting two plaquettes: e.g. 1-1 and 2-2 in Fig. m 
and compare it to the pair of hops 1-0 and 0-1 in Fig. [31 
(Indeed, to make the tube into a ladder, we could cut the 
inter-plaquette tube bonds 0-0 and 3-3, alternating with 
cutting 1-1 and 2-2, down the line, then flattening it out 
like an accordion-fold.) Most importantly, looking in par- 
allel directions along the adjoining edges, the numbering 
around the B plaquette has the opposite sense from that 
of Fig. m We chose to use the numbering scheme shown 
so as to keep the clock sense the same on both plaquettes. 

Thus common sense combined with the pseudospin 
dictionary (j2.5p explains why the x and z components 
flip sign, but the y component does not, in p.7p . For 
example, imagine an interplaquette interaction that fa- 
vors having fermions at both ends of a weak bond (see 
Sec. IIIIDI below). We see in that for the tube (Fig. ^ 
that means (say) fermions take even sites on both pla- 
quettes, but the same thing on the ladder (Fig. [3]) means 
they are even on one plaquette and odd on the other; the 
difference between even and odd charge-order is a sign 
fiip of the Px or component. Also, the sign of pseu- 
docurrents (and hence of Pz or Tz) is manifestly flipped 
when the sense around the plaquette is reversed. On the 
other hand, if an interaction means that (say) bond-order 
on the 1-2 bond of one plaquette of the tube (in Fig. [3]) 
repels bond order on the 1-2 bond facing it, i.e. favored 
antiferromagnetic alignment of the Py or Ty components, 
the same thing is true for the 0-1 bonds on the ladder 
plaquettes (in Fig. 2]) . 

In Eq. (|3.8p . there is a uniform magnetic field in the 
pseudospin y direction. We get the single-site terms 
from the "excited" hops (defined at the start of this sec- 
tion) which are not disallowed for the ladder. This com- 
petes with the antiferromagnetic exchange term, having 
the effect (as usual in antiferromagnets) of a uniaxial 
anisotropy favoring the xz plane. Hence, the system has 
the symmetry of an XY model ordering in that plane. 
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which corresponds [by Eq. (|2.5p ] to CDW and sponta- 
neous currents. Similar to the tube, having a continuous 
symmetry in one dimension, it would only have power- 
law correlations. 

Doing the same for the square lattice amounts to ex- 
tending the result of the ladder calculation to a square 
lattice. Recalling for the ladder [Eq. (|3.8p ]. the pseu- 
dospin Hamiltonian for the two-plaquettc unit was 



= —[const{x) +Ta-Tb]^ — ' 



V 



2V\l + x 



(3.9) 

The similar result for the perpendicular direction in the 
plane would be 



t'- 



■[const{x) -\- Ta ■ Tb] 



t'- 



.(- 



1 



(Tl + Tl) 



(3.10) 

The minus sign for the single-plaquette terms in the sec- 
ond case is because the bond-ordering in the two per- 
pendicular directions are the pseudospin in -fy and ~y 
directions respectively (see FigH]). Hence for the infinite 
square lattice, we get 



12 



Hsquare = 2^ [cOnst{x) + Ta ■ Tp] 



(3.11) 



which is the antiferromagnetic Heisenberg Hamiltonian. 
Since, the square lattice is two-dimensional, its ground 
state will possess long range order. Notice that antifer- 
romagnetic tendency of staggered pseudospin implies a 
ferromagnetic tendency for spontaneous currents in both 
ladder and the square lattice. 



C. Role of Symmetries 

Exploring symmetries can lead to a better understand- 
ing of the relation between the form of microscopic model 
and that of the effective Hamiltonian. Our starting 
fermion Hamiltonian (acting on a a two-plaquette unit) 
had the following symmetries: a) time reversal b) reflec- 
tion symmetry (flipping the two plaquette upside down). 
We have only considered models that maintain these 
symmetries. 



1. Consequences of generic symmetries 

These symmetries imply specific symmetries in the 
pseudospin effective Hamiltonian (|3.ip . (a) The absence 
of a single-site term follows from the microscopic 
time-reversal symmetry, under which the the pseudocur- 
rent operator flips sign. (b). The absence of a single-site 
term is due to the transverse reflection symmetry of 
the two-plaquette unit, under which the charge-density 
operator flips sign.^^ 

The tube and square lattices both have a 4-fold ro- 
tation symmetry, too, under which the two bond-order 



states (pseudospin +y and —y) are equivalent, ergo the 
py terms are absent. Morever, for the most general 
one-particle spectrum that a single plaquette could have 
(keeping intact the degeneracy of momentum carrying 
states, i.e., -1-1 and —1, but lacking particle-hole sym- 
metry) , the effective Hamiltonian would still be of form 
Eq. (j3.8p . though with different numerical coefficients. 
On the other hand, the effective Hamiltonian for the lad- 
der generically includes single-site P^ terms, since they 
are not ruled out by any symmetry. 



2. Role of lattice symmetry in ladder model 

As we just discussed, a four-fold lattice symmetry guar- 
antees certain pseudospin symmetries. In our basic lad- 
der Hamiltonian (|2.ip , the single-plaquette terms had an 
"accidental" (non generic) four-fold rotational symmetry 
not guaranteed by the ladder's symmetries. If we general- 
ize the ladder model so as to break the fourfold symmetry, 
what kinds of pseudospin asymmetries are generated? 




FIG. 5: A ladder without four-fold rotational symmetry 
around a plaquette. 

First we can make the transverse hopping tj_ within a 
plaquette different from the longitudinal t (Fig. [5|) . Then 
the effective Hamiltonian turns out to be 



t'- 



Hladder =COnst + X! ( OT? ) ' "^"+ 



2V 
2At 



2V J \ {V -2t + At)iy -2t + 3At) 



rpyrpy 



2(y-2i + 3Ai) ^^^V^" 



(3.12) 



where t = (< + tj^)/2 and At = — t. The absence of 
four-fold symmetry of the hopping around the plaquette 
leads to a first order field term -|-2AiT*' and a second- 
order anisotropic exchange along the Ty (BDW) pseu- 
dospin direction for the ladder. The underlying reason 
for this is that, so long as = ij^, angular momentum 
is a good quantum number and in the m — ±1 single- 
fermion states, the sites from which fermions can hop 
differ in phase by ±7r/2; the upshot is that angular mo- 
mentum is conserved by degenerate hops, even though 
(in the ladder) it is not conserved by excited-state hops, 
and thus the exchange terms are isotropic. 

Rather than spoil the ladder plaquette's four-fold sym- 
metry in the hopping terms, which couple to the pseu- 
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dospin Ty component, we could do it by making intra- 
plaquette interactions unequal, which produce single- 
pseudospin terms. Making transverse and longitudi- 
nal interaction different gives a linear coupling at first- 
order in perturbation theory to the bond-charge operator 
(X (V^i — V±)Ty; as before, the time reversal and reflection 
symmetries forbid linear or terms >^ 

D. Engineering spontaneous currents by fermion 
interactions (ladder) 

Our study was motivated by the question : Can the 
microscopic models considered so far exhibit current car- 
rying states spontaneously in their ground state with gen- 
uine order in an Ising sense. In the language of the pseu- 
dospin mapping [Eq. (|2.5p . we want the effective Hamil- 
tonian Eq. (j3.ip to give spin order along a particular 
axis. But as we saw in Subsecs. IIII Al and IIII Bl the 
natural form of (|3.ip had a continuous XY symmetry 
(corresponding to the currents/charge operators) in the 
ladder case; in the case of the tube or square lattice, 
(j3.1|) had the full three-component rotational symmetry 
of a Heisenberg magnet. To stabilize any particular kind 
of order more than the others, we must spoil these un- 
wanted symmetries by anisotropic pseudospin terms hav- 
ing the effect of an Ising-like anisotropy in the desired 
direction. In this subsection and the next, we continue 
the preceding one by surveying various generalizations in 
the fermion Hamiltonian and and the pseudospin terms 
they give rise to, but now motivated by a sort of engi- 
neering: rather than just solve for the ground state of 
a given Hamiltonian, we frankly seek the Hamiltonian 
which gives the targeted state. 

That need not require getting a uniaxial easy-z 
anisotropy (which - see Subsec. IIII E II and IIII E 21 - is 
impossible from the fermion Hamiltonians we admit). 
Instead, we look for two perturbations that combine to 
disfavor the other two directions. Now, since we have 
(pseudo)spin-l/2, there is no way to get any uniaxial 
anisotropy from a single-spin term, but a trick is avail- 
able: given we have antiferromagnetic order, a uniform 
external field perturbation creates an effective easy-plane 
for the directions normal to the field; this was already 
used to get xz anisotropy for the ladder case (from the 
linear term in (j3.8p V The other perturbation giving 
the needed hard-x uniaxial anisotropy must take the form 
of a pseudospin-pseudospin interaction, coming from the 
inter-plaquette terms of the fermion Hamiltonian. 

One natural extension of our model is to add an inter- 
plaquette nearest-neighbor interaction 

H' = V'{nA{l)nBiO) + ?iA(0)ns(l)) (3.13) 

We have used the ladder numbering scheme in writing 
the above expression. Looking at Fig[2l we see how the 
term p.l3p distinguishes the CDW sector from others, 
since the operators like 71^(1) depend on the CDW order; 
on the other hand, it cannot distinguish different CCS 



or BDW states, since they have equal fermion densities 
on all sites. Using Eq. (|2.5p . we can easily convert the 
interaction term to pseudospin language and indeed 

, /I PaPr\ //I TaT^\ 
K -V'[-,- ^) =V'[-, + ^) (3.14) 

We emphasize the above effective interaction p.l4p is 
first order in perturbation theory, and not second order 
as for the hopping processes earlier in Sec. lIIII Eq. (|3.14p 
is an adjustment of the Jx pseudospin coupling in (j3.ip 
and thus favors CDW order, either a uniform pattern 
on each plaquette or an alternating one, depending on 
the sign and magnitude of V'. An inter-plaqutte sec- 
ond nearest- neighbor interaction , i.e. y'(n^(0)nB(0) -I- 
nyi(l)ns(l)), gives the same result as (|3.14l) but with a 
flipped sign for the exchange term. 

Thus, we see a route to favoring spontaneous currents 
for the ladder. An infinitesimal attractive inter-plaqutte 
second nearest-neighbor interaction (V' < 0) or repul- 
sive inter-plaquette second nearest-neighbor interaction 
{V" > 0) will make Jz > Jx and the ground state will 
have currents spontaneously. For a comparison, we note 
that in generalized Hubbard models, attractive nearest 
neighbor interaction was argued to stabilize currents^^. 

For the tube and square lattice, inter-plaquette inter- 
actions can only reduce the antiferromagnetic Heisen- 
berg symmetry to a continuous XY symmetry in cur- 
rents/BDW plane and do not favor currents exclusively. 
To do that, we must look to interplaquette hoppings in- 
stead. 



E. Engineering spontaneous currents by fermion 
hops (square lattice) 

An alternative extension of our model is to add addi- 
tional inter-plaquette hoppings. As we will see, for the 
ladder and square lattice, this favors bond (BDW) order 
by increasing the Jy pseudospin exchange while decreas- 
ing Jx and Jz couplings in (|3.ip . For the tube, inter- 
plaquette hopping to any distance can never reduce the 
continuous Heisenberg symmetry due to 4-fold symme- 
try. 

1. Guessing the fermion term? 

A short-cut may allow us to quickly find fermion terms 
yielding a desired inter-plaquette pseudospin Hamilto- 
nian form. Let's extend the notion of "pseudospin" back- 
wards to impute pseudospin to the single-fermion states 
I ± 1). Indeed, we can just ignore the other orbitals, 
since only the "degenerate" hoppings (explained at start 
of Sec. imp could give us a pseudospin interaction from 
second-order perturbation theory. Then, we just substi- 
tute Ta — > J2a'a '^lr'i^)^'^'<^'^<^i^)' whcrc TA mcans the 
usual Pauli matrices. Thus, any coupling Ti-^^ff between 
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(components of) pseudospins Ta and Tb gets transcribed 
to a four-fermion term Ti.^g. If we can regroup the four 
fermion operators so that H^g cx —H'^^H' where the op- 
erator H' hops a fermion from plaquette B to A, then we 
could take Ti.' to be the inter-plaquette term reducing to 
Ti.'^g via second-order perturbation theory. 

When apphed to two plaqucttcs in a ladder or square 
lattice (using the site labels of Fig. 2]), we get 

Ta^b iCA,0 - c\ qCa,i) ■ i{cB iCBfi - Cb oCb,i) 

' (3.15) 

and (for the bond-order component) 

(3.16) 

Thus, 

nn+nn ^ (c^,iCB,o)(cA,ocki) +{a^b) (s.n) 

where there was only one of grouping such that the inter- 
plaquette hops connected nearest neighbors; this is the 
term we found already. What about the term which 
would break the degeneracy between bond order and 
pseudocurrents : 

- Txn + nn - (cLocb,o)(ca,i4,i) + {a^b) 

(3.18) 

The grouped factors in Eq. (|3.18p are diagonal ("iy2") 
fermion hops. Unfortunately, the sign of this term is 
necessarily positive, so it always favors the Ty (bond- 
order) direction. 

2. Fermion hops in general 

A more comprehensive study of hops will be prof- 
itable for the following reasons: (i) As we are about to 
show, it reveals that the findings in (|3.17p and (|3.18[) 
are general for any nearest-neighbor hopping, so that 
is not a route to the desired order, (ii) Consequently, 
further-neighbor pscudospin interactions coming from 
long-distance fermion hops are our last hope to disfavor 
bond-order in the square lattice, and the general formula 
guides us to the correct interactions for this purpose (iii) 
It is the root reason that the ladder and square lattices' 
pseudospins needed to be staggered (at the beginning of 
Sec. imp but not the tube lattice's. 

So imagine a perturbation Hamiltonian containing 
hops from any vertex of one strong plaquette to any 
vertex of another (not necessarily the nearest neighbor). 
First consider (still) the nearest-neighbor plaquettes, and 
let t'l^ i'yj, t'2, etc. be (weak) hoppings to sites (in the 

other plaquette) at distances 1, v^, 2 and so on, re- 
spectively. Then the exchange part of the two-plaquette 
effective pseudospin hamiltonian is 

= JtTa ■ Tb + JpPa ■ Pb 

= {Jt + Jp)PlPl {Jt - Jp){P1Pb + PIPb) 

(3.19) 



where 

Jt - (t'l - 2t'^ + t'^f/V > 0, (3.20a) 
Jp = {t'^ - 2t', + t'^f/V > (3.20b) 

Which term does a given hopping contribute to? First, 
remember we always need two different hoppings, with 
endpoints distinct; if they end at the same site on one 
plaquette, we could not couple to that plaquette's angu- 
lar moment. Then if we orient the two plaquettes such 
that the two hoppings don't cross, the exchange coupling 
relates staggered pseudospins or plain pseudospins de- 
pending on whether clock sense on the two plaquettes 
arc same or alternating respectively. E.g., in Fig. |4l 
the non-crossed hoppings are connecting plaquettes with 
same clock sense, thus giving rise to staggered pseudospin 
exchange coupling, i.e. Ta ■ Tb- 

A related observation is that, in the ladder, there is a 
symmetry under mirror-flipping every second plaquette 
around the long axis of the ladder, while switching t[ <-!■ 
^V2 ^^'^ ^ ^'\/5' (Ei^ni) shows this switches the T 
and P terms. 

Inspecting (|3.19p . we see that so long as we have only 
crossed or only uncrossed fermion hoppings, the result is 
isotropic in the {yz) pseudospin plane, so that bond order 
and currents are degenerate. However, if we start from 
a mixture of crossed and uncrossed hoppings, the bond- 
order (y) exchange is always stronger than the currents 
(z) exchange — and is always antiferromagnetic. 



3. Spontaneous currents via anisotropic frustration 

Given this last fact, is it possible at all to obtain a pseu- 
dospin anisotropy favoring T^ and hence current order 
over the whole lattice, by coupling more distant units? 
This is possible, in principle, through anisotropic frustra- 
tion. (It is assumed interactions have somehow already 
disfavored charge ordering, as discussed in subsubsection 

nmi) 

Assume the dominant nearest-neighbor hopping is 
purely t'l, as in our original and simplest model. The 
pseudospin exchange has a continuous symmetry in the 
yz spin plane, leading to antiferromagnetic order degen- 
erately in any mixture of those components. Now imag- 
ine (say) a second-neighbor pseudospin exchange due to 
mixed kinds of hoppings; by the above arguments, J2y is 
necessarily antiferromagnetic, and J2y > \ J2z\- But un- 
like the nearest-neighbor exchange, the enhanced second- 
neighbor J2y term disfavors bond-order state (being of 
the wrong sign). 

A second-neighbor exchange is allowed on the ladder, 
using (say) the hops t"^ and t" connecting two pla- 
quettes related by a [4, 0] vector. On the square lattice, 
the second nearest neighbor has a displacement [2,2] and 
this exchange turns out to be symmetry-forbidden. But 
the square lattice can have the same [4, 0] inter-plaquette 
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hops as on the ladder, and these finally give our goal: we 
can favor spontaneous order in the square lattice, albeit 
with a rather baroque Hamiltonian. 



IV. GENERALIZATION TO SPINFULL 
MODELS 

It is natural to ask if we can extend our results to 
models with spin. In case they can be applied to a real 
electronic system, and also to make some contact with 
other works on spontaneous current models. There are 
two quite different ways to imagine this. First, as worked 
out in Subsec flV A| we can simply include an additional 
spin degree of freedom in the Hamiltonians considered 
above. Alternately, as worked out in Subscc. IIVBI we 
can exactly map a site degree of freedom in one of our 
spinless models to the spin degree of freedom in a model 
with half as many of lattice sites (thus keeping constant 
the the total degrees of freedom.) 



A. Adding Spin Degree of Freedom 

For this extension of our model, we simply add spin 
indices in all the terms of (|2.ip while conserving the spin, 
and rerun the calculations of Sees. [TTl and IIIII 

To make the spinfuU calculation analogous to what we 
did, the filling should now be 3/8. Of the three fermions 
per plaquette, the first two fill angular momentum zero 
(|0 t) and |0 |)). The third fermion goes in the de- 
generate current carrying state, | ± la). The extended 
pseudospin representing each plaquette is now the direct 
product of the same pseudospin degree of freedom, and a 
real spin Sa —■ Also, the interaction term in Eq. (|2.1b[) 
has to be augmented to include an onsite interaction term 
([/) equal in strength to the offsite interaction terms {V) 
so that the multi-fermion eigenstates still remain direct 
products of single- fermion orbitals. 

Here are the results for each case: 
SpinfuU Tube : 



+4(F„-F^+i)(5„-^„+i)l (4.1) 



SpinfuU Ladder : 



Heff - 



1 + 



(4.2) 



SpinfuU Square Lattice : 

^^ff = (^)e [Q)^-^^ 



/3 ■ 



{Ta ■ Ti3){Sa ■ Sp) 



(4.3) 



{Ta ■ Ta+l){Sa ' Sa + l) — 



1+X 



Thus the effective Hamiltonians have a form like the 
Kugel-Khomskii Hamiltonian^l for cubic titanates de- 
scribing spin and orbital superexchange interactions be- 
tween ions having threefold degenerate t2g orbitals. 

The result (|4.ip for the tube case is proportional (mod- 
ulo a constant) to {^+2Pa-Pa+i) {^+'2Sa-Sa+i) which is 
the SU (4) symmetric Kugel-Khomskii modeP^ . For the 
tube, the interaction terms are just "degenerate" hops of 
Sec. IIIII they conserve spin as well as pseudospin. Actu- 
ally, they conserve a combined fiavor which includes both 
the spin and the pseudospin. Hence, the effective Hamil- 
tonian possesses an SU (4) symmetry in which there is no 
distinction between the four combined flavors the hop- 
ping fermion might carry. The interes ting behavior of 
such SU (4) chains is discussed in Ref. [33; in terms of 
the original fermions, it obviously corresponds to a high 
degeneracy between many kinds of order. 

For the ladder and square lattice cases, the "degener- 
ate" hops do not conserve the combined flavor thereby 
reducing the S'C/(4) symmetry to only SU{2) x SU{2) for 
the exchange terms. What kind of order do these lat- 
tices have? Notice that the spin-pseudospin cross-terms 
tend to favor ferromagnetic order in one kind and anti- 
ferromagnetic order in the other. Since we also have a 
pseudospin antiferromagnetic exchange but no real spin 
exchange, the expected order is always ferromagnetic for 
the real spins^. Then the pseudospin order is the same 
as in a spinless model; in effect, the system spontaneously 
becomes spinless by polarizing in one spin flavor. 

We compare our results to that of'^^, in which a Hub- 
bard model with a similar pattern of strong and weak pla- 
quettes {t and t') with just an onsite interaction term was 
studied on a square lattice. Yao et al found a host of dif- 
ferent phases for different values of the onsite interaction 
including a Fermi Liquid, a d-BEC, a d-CDW, a d-BCS, 
a spin-1/2 antiferromagnet, a spin-3/2 antiferromagnet 
and an "orbital nematic" phase at 3/8 filling, which is 
equivalent to one doped hole, with respect to half filling, 
on each plaquette ("Q^t = 1" or "a; = 1/4" in their nota- 
tion). In their model, it is only in the parameter regime 
Uc — 4.6i < U < Ut — 18. 6i that their single-plaquette 
states are filled like ours and admit the possibility for 
currentsSS Then, the single-plaquette states are charac- 
terised by spin-1/2 as well as a pseudospin-1/2 (caUed 
chirality Tz by Ref. [H, and having "p^ ± ipy" symmetry 
i.e. our angular momentum | ± 1).) 

They do not discuss the possibilities of orbital currents 
explicity in the above mentioned regime, but implicit in 
their result are the anisotropics in the pseudospin ex- 
change terms which is interesting (See Eq. (4) and (6) 
oS§.). They state that, in the regime Uc <U <Un — 7.3, 
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the system becomes a spin- 1/2 antiferromagnet with 
electron nematic order (same as our BDW), while in 
Un < U < Ut, there is no nematic ordering. Perhaps, 
there are spontaneous currents in this regime. How- 
ever, they do not discuss the origin of the pseudospin 
anisotropics. It is all the more perplexing to us, given 
our experience that one set of "degenerate" hops taking 
part in the second-order perturbation theory can only 
give rise to isotropic exchange. 

B. The Nersesyan Map 

As first proposed by NersesyaniS, a spinfuU model on 
any lattice can be mapped to a spinless model on a dou- 
bled version of that lattice (its direct product by {1, 2}.) 
Each pair of sites in the spinless model represents re- 
spectively the spin-up and spin-down occupation. Thus, 
a spinless ladder maps to a spinfuU chain (or vice versa), 
such that the leg index maps to the spin index. (We shall 
call this a "rung spin" to distinguish it from the real spin 
of Sec. IIV Al and the pseudospin of all the earlier sec- 
tions.) Hamiltonian terms acting on a single rung of the 
ladder will be mapped to single-site terms on the chain, 
while terms along the ladder's leg map to terms along 
the chain's leg. We exhibit examples of the map in both 
directions. 

In fact, since our plaquette is built from two rungs, 
each plaquette pseudospin operator Pi corresponds to 
two neighboring rung spin operators as shown below. 
Consider the fermion basis states on each rung j that 
have nonzero pseudocurrents, namely 

10,+) = (|j,l)+*|j,2>)/\/2, (4.4a) 

^ (|j,l)-*b-,2))/V2. (4.4b) 

The ± label in (|4.4|) is a rung pseudospin index defining 
the z-axes for rung pseudospin S' j aligned with that of 
the plaquette pseudospin such that 

- {P^I^.P!) (4.5a) 
S'2j+i = (-P/,P^j,-P/) (4.5b) 

On the other hand, Nersesyan's rung spins (we call S 
keeping in mind the difference from the previous section) 
are related to our rung spins S' via (See Eq. 3-5 o^) 

= [-lys'j' (4.6a) 
= {-lySf (4.6b) 
5| = [-lySf (4.6c) 

Notice that the definition of Nersesyan's rung spins is 
staggered compared to ours. Nersesyan used a "canted" 
rung spin basis to make manifest the staggered na- 
ture of charge/spin densities in CDWs and Spin Den- 
sity Waves(SDW) for a spin-1/2 Hubbard chain. We will 



write rung spins using the Nersesyan basis from now on. 
Using (|4.6p . we can re-use the results of Sec. HIT] to fig- 
ure out what a given fermion perturbation projects to in 
terms of Nersesyan pseudospins. 

We shall return to apply the Nersesyan map in 
Sec.lVB] 

1. Spinless ladder to spinfuU chain 

In particular, our basic Hamiltonian for the spinless 
ladder [see Eq. (|2.ip ] gets mapped to a chain of alternat- 
ing strong and weak bonds ("sf" here distinguishes the 
spinfuU model parameters): 

n = - Y,tf^,+ MaC^+l■- + h.c] 

i i i 

with /if = t, C/"f = V; we get = t or t', and 

Vfl^i = or 0, respectively for strong or weak bonds 
+ Thus our spinfuU chain includes the usual hop- 
ping and interaction a nearest neighbor interaction and 
a field term along x in spin-space or a single-site spin flip 
term. This kind of spin mapping will map spontaneous 
current states to spin current states with canted site-spin 
expectations, CDW to SDW, and BDW to an equivalent 
BDW/paramagnet as one may readily verify. The advan- 
tage of this kind of mapping is that we may carry known 
results from spinfuU models to spinless models or vice 
versa. 



2. SpinfuU ladder to spinless tube 

Tsuchiizu and Furusaki"^- presented a spinfuU ladder 
model which had a spontaneous current phase. Their 
Hamiltonian included the following kinds of parameters 
Two hoppings (longitudinal i^^^y and rung t^^±); an on- 
site Hubbard repulsion U^^; three neighbor repulsions 
(ladder V^^^ , rung V]_^ , and second-neig hbor F^f); a spin 

exchange J+, acting across a rung; and, last but not 
least, a correlated hopping or ring exchange t^^^j., which 
takes two electrons from one site to the other site on the 
same rung. 

Under the "Nersesyan map" , their model corresponds 
to a kind of "tube" model with spinless fermions, but of 
course lacking the 4-fold symmetry of our tube model. 
Their repulsions V^^ and U^^ , as well as the zz term of 
the exchange interaction, map to various on-plaquette 
repulsions like out V; while their Vj^* and V^^ map to 
combinations of inter-plaquette interactions in the spin- 
less model, including (but not limited to) our V' . Their 
tf is the same as our inter-plaquette hopping t' , but 
their rung hopping t+ corresponds to hopping on only 
two edges of the plaquette in the spinless model. (An 
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x-oriented magnetic field - not included in their model - 
would have mapped to spinless hopping on the plaque- 
tte's other edges.) Finally, their spin-exchange and 
correlated hop t^^^^. both map to correlated (two-fermion) 
hops, like no term in our models: the xy component 
terms of would correspond to two diagonal spinless 
fermions on one plaquette moving to the opposite diag- 
onal, while ip^j^jj. takes two adjacent fermions on one side 
of a plaquette and moves them to the two opposite sites. 

The strong coupling approach of Refi^ corresponds 
exactly to our separation of "weak" and "strong" pla- 
quettes. Their state with spontaneous (pseudo)currents 
is somewhat disappointing, as the return current depends 
on their pair hopping term. (Of course, in this limit that 
all inter-rung couplings go to zero, it would not be pos- 
sible to have a genuine current involving two sites, since 
conservation requires a loop for it to circulate around.) 

SchoUwock et aW^ exhibited a simpler spinfuU ladder 
which was shown numerically to have a spontaneous- 
current phase. They have just one hopping parameter 
f'l = = t^^ (thus they could not access the strong- 
coupling limit of small t||^). They also have Hubbard 

U , exchange Jj^, and repulsion Vf; . In the special case 
that J^^ 1^= U^^ — V]^, this maps to a tube model with 
fourfold symmetry (this is not the same special point 
Jst/4 ^ [/St + vf, which had S'0(5) symmetry^). The 
mapped spinless model would also have the same corre- 
lated hop of diagonal pairs as Ref. [H. 



V. TOWARDS UNIFORM SPINLESS MODELS 

Our hope was that if we can find a strong-coupling 
model that has spontaneous currents, perhaps it can 
be adiabatically continued to a translationally invariant 
model that does not distinguish the Hamiltonian terms 
on weak and strong plaquettes. There are two precondi- 
tions to even consider this: 

(i) The symmetry pattern of order should be con- 
sistent with a uniform order. Below, in subsec. fV Al 
we verify this for our models. 

(ii) The "weak" and "strong" terms in the Hamil- 
tonian should have the same form, differing only 
by the size of the coefficient. That is easy enough 
to manage, even if we adhere to the somewhat un- 
natural interaction term (|2.1b[) . with the nearest- 
and second-neighbor strengths made equal for con- 
venience. If that equality is carried over to the 
inter-plaquette terms, it is actually beneficial since 
it cancels a term favoring CDW ordering [see Eq. 

(EH] 

There is one further challenge: having conjectured a 
Hamiltonian favoring current ordering, how could we ver- 
ify that? We need some family of variational wavefunc- 
tions that would (ideally) be definable for all the interpo- 
lating Hamiltonians from strong coupling to uniform, and 



where the variational parameters allow any value of the 
spontaneous-current order parameter. (The best model 
of such a calculation is the recent work of Ref. H.) In the 
case of the spinless ladder models, such a correspondence 
was already worked out by Nersesyarti^ and is elaborated 
below in Sec. IV Bl 



A. Inter-plaquette Spontaneous currents 

Previously, we saw that at least for the ladder, the 
ground state could exhibit spontaneous currents by pro- 
ducing anisotropy along that sector. Such a symmetry- 
broken state has spontaneous current expectations not 
only on the strong plaquettes, but also on the weak bonds 
connecting them. Depending on the symmetry of this 
pattern, it may be possible to adiabatically connect to 
a state in which the strong/weak distinction goes away 
and the current expectations have equal strength on all 
bonds. 

Current expectations on Weak Bonds: For the 

weak bonds, the required pseudo-current operator is 
iweak = i[c\iCBj - c^B^'Ai] whcrc i - j IS weak bond. 
We again use the canonical transformation recipe, (|A2[) 
from the appendix, but now for the purpose of deriv- 
ing an expectation rather than a Hamiltonian term. The 
pseudo-current operator for the weak bonds projects, in 
pseudospin language, to 



^weak 



t' 



2(2t 4 


-V) 




2-) 


t' 




2(2t-f 


-V) 


t' 




2(2t 4 


-V) 



(|2+;2+)(2+;2- 



{Pl + Ph) 



in - n) 



(5.1) 



Thus, in the ladder or square lattice where Pz tends to be 
uniform, i.e. is staggered, we see the (pseudo)currents 
in weak bonds are proportional to those in the adjoining 
strong bonds: both kinds of current must have the same 
degree of order. However, in the tube lattice, the weak 
bond currents are zero. That must be true to all orders 
on grounds of symmetry. (Reflection in the plane of a 
strong plaquette reverses the sense of currents on weak 
plaquettes but not on strong ones.) 

Since the pseudospin operator for the strong bond is 
just /2 [Eq. 12.5] . the weak-bond and strong-bond 
currents have opposite directions (Fig. [5]), which is the 
same pattern as the c?-density wave stated. Thus, this 
pattern is consistent with a ddw-like state if we could 
analytically continue our model to a uniform one. 



B. Uniform model 

This is not the same as our pseudospin map of Sec. Ill Bl 
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FIG. 6: The orbital current pattern in |2+; 2+) state. 



Can one of our strong-coupling models be connected 
to a uniform one, in which the distinction of strong 
and weak bond vanishes? We found the Nersesyan map 
[Sec. IIV B| was the key to guessing a Hamiltonian that 
has a spontaneous-current ground state. The spinless 
model's pseudocurrent order maps (as we show shortly) 
to a certain spin order in the spinful model. One need 
only invent a uniform spin Hamiltonian giving the de- 
sired order, and such that it maps to a plausible Hamil- 
tonian of the spinless fermions. The only limitation is 
that the resulting spinless model is built from site pairs: 
it lends itself to ladders, or bilayers, but not to (say) the 
fourfold symmetric square lattice. 

The Nersesyan approach is an existence proof that if we 
get spontaneous current order in a simple, local, strong- 
coupling model, the adiabatic extension can work all the 
way to the uniform case. If we could find an analogous 
starting point on the square lattice, that result would 
make the corresponding extension more plausible; unfor- 
tunately (Sec. IV B ip we could not find such a starting 
point. 

Concerning the certain spin order mentioned above, 
consider a strong plaquette in our ladder with angular 
momentum -f 1 which has spontaneous pseudocurrents on 
its rungs and leg bonds. It maps to a pair of sites on the 
spinfull chain with Nersesyan spin expectations 



(5,)(x(l,(-l)^o) 



(5.2) 



which corresponds correctly to a non-zero (PJ) — 1 (See 
Eq. g31) and Eq. gS])). The unifo rm Sx component 
is not surprising, as every spinless ladder we consider 
includes a rung hopping which becomes an x transverse 
field under Nersesyan's map. 

The pseudocurrent along a rung, I^, maps to Sy, cor- 
responding to the alternating y component. What about 
the pseudocurrents along the legs? Let's define a differ- 
ence between the two sides, 



I (i, 1)^(1+1.1) - 2)^(1+1,2)- 



(5.3) 



Of course, in the expected state, opposite sides have op- 
posite currents, so this also has a nonzero expectation. 
Evidently, this simply maps to the z component of spin 
current from i — > i + 1, Indeed, in spin mod- 

els with an isotropic Hamiltonian (and also here), spin 
current goes with the twist between noncoUinear spins, 
li^j oc Si X Sj and indeed, the staggering of rung-spin 
directions in (|5.2p does give the requisite nonzero (and 
alternating) (Si x Si+i)z component. 



Now we see a simple route to a uniform model having 
spontaneous currents: simply find a uniform Hamiltonian 
for a spin- 1/2 chain, having a ground state with the above 
staggered spin canting. The simplest realization (from 
the spin chain viewpoint) would be an antifcrromagnetic 
chain with (i) an Sx field, so the spins will cant transverse 
to it - recall this is a simple transverse hopping in the 
spinless model; plus (ii) a small anistropy in the antifcr- 
romagnetic exchange, such that the Sy axis is easier than 
Sz] this ensures the canting happens in the y direction. 
The problem with following this route literally is ingre- 
dient (ii): Sx-Sx or Sy-Sy spin couplings correspond, in 
the spinless model, to two-fermion pair correlated hops 
in a plaquette, which we did not want to include. Con- 
versely, this nicely illustrates why correlated hopping is 
conducive to the existence of spontaneous currents^ ! 

The scenario of the previous paragraph can be 
achieved, instead, with the following adjustment of the 
Hamiltonian: in place of a spin chain, take a Hubbard 
chain with a transverse x field h^J , as above, plus a 
small ferromagnetic coupling of neighboring S^ com- 
ponents. A strong Hubbard on-site repulsion U^^ pro- 
vides the effective antiferromagnet exchange at order 
[t^^Y in the standard fashion. Via the Nersesyan 
map (see Eq. (|4.5p and (|4.6p ). we realise that: 1) the 
transverse field suppressing spin ordering in x direction 
maps to a term suppressing of BDW in the spinless lad- 
der, 2) the small ferromagnetic z coupling suppressing 
spin ordering in z direction (by reducing the strength of 
z antifcrromagnetic exchange) maps to a term suppress- 
ing CDW in the spinless ladder, and finally 3) due to the 
aforesaid suppresions, stabilization of spin order in the 
y component maps to orbital currents for the spinless 
ladder. 

In the spinless language, the Hubbard U^^ and 
—jfSfSfj^i interactions just map, respectively, to 
nearest-neighbor interactions along the rungs and legs 
(repulsive and attractive V||). The Hubbard hopping 
t^^ and the transverse field just map, respectively, to 
hoppings along the rungs and legs (our t±^ and <||). The 
half filling we adopted for the spinless model corresponds 
to half filling in the Hubbard model. Unlike all spinless 
models we previously mentioned, this model is a uniform 
ladder with no weak and strong plaquettes. 

The above paragraph is essentially a rediscovery of 
Nersesyan's ladder model, Eq. (1) of^-. In particular, 
our key ingredient - making V|| attractive while V± repul- 
sive and zero - is essentially the same as Nersesyan's 
recipe, which is that K/2— V|| > 0. It is interesting to note 
these chains have fractionalized excitations, domain walls 
carrying 1/2 fermion charge^S, corresponding to spinous 
in the spin model. 



1. Extension to the square lattice? 

The same mapping can be used in one dimension 
higher, to build a spinless spontaneous-current state on 
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a bilayer from a Hubbard model on a square lattice^i 
(or any bipartite lattice). A spinfull model with plau- 
sible interactions on that same lattice was known ear- 
lier— that exhibits spontaneous currents. That model's 
Hamiltonian is similar to the SchoUwock et al ladder"^^ 
we described above, in one dimension higher, and is sim- 
ilarly engineered to have an S0(5) symmetric point. (A 
minor difference is their "rung" (interlayer) hopping 
may differ from the in-layer hopping t||^.) All their inter- 
action terms - , Vf, and (isotropic) exchange jf^ - 
act only on rungs. Thus, the main qualitative difference 
our model differs (apart from spin, of course) is our in- 
layer repulsion term Vjj . The (conjectured) order in our 
bilayer model is the same alternating pattern of currents 
as in their model (Fig. 1 of Ref. issh . 

Unfortunately, this does not work for the plain square 
lattice, for two reasons^. Firstly, the best we could man- 
age for a Hamiltonian fSec. IHI E 3p depended on t" and 
^"yio hoppings to a third neighbor plaquette, while the 
much shorter hoppings entering (j3.20b[) must be negli- 
gible. Not only are those absurdly distant neighbors to 
have a meaningful hopping: once we make the lattice 
uniform, we must also include t' ^ and t'^jo same 
separations as t" ^ and t" but coupling nearest neigh- 
bor plaquettes and creating the undesired un-frustrated 
anisotropy (|3.20bp . Secondly, we cannot be guided by 
Nersesyan's map, as it demands that sites come in pairs. 

Although, it seems difficult to obtain currents order 
in a uniform square latttice, the terms which worked in 
the ladder would still be effective in a rectangular lattice. 
That suggests one possibility to obtain spontaneous cur- 
rents as a symmetry breaking. Notice that a uniform 
bond order, i.e. the state currently called "electron ne- 
matic"— reduces the system to rectangular symmetry, 
making t\\ ^ t_\_ heuristically (See Sec. IIIIC2[) . Thus, 
spontaneous currents could be parasitic on electron ne- 
matic order. This is a not a linear coupling of the two 
order parameters; it would be a second Ising-like transi- 
tion, to a state breaking Z2 x Z2 symmetry. 

VI. DISCUSSION 

The central aim of this paper was to investigate the 
possibilities of a lattice model manifestly displaying spon- 
taneous currents in its ground state. Among our mod- 
els, containing standard hopping and interaction terms, 
it was fairly difficult to stabilize only the spontaneous- 
current state. Frequently, there was a remnant contin- 
uous symmetry leading to an arbitrary mixing with one 
of the competing orders; and most perturbations which 
could break that degeneracy tended to favor the com- 
peting order. Something similar also happens in some 
spinfull models meant to address the possibility of spon- 
taneous currents in a realistic system, e.g. the relation2£ 
of "d-density wave" current order and "electron-nematic" 
order (related to our "bond order" ) . 

We did show that rather contrived and unappealing 



fermion interaction or hopping terms could be used to 
stabilize currents, but it seems highly unlikely that those 
kinds of processes are at work in real materials. We sug- 
gest our finding may be related to the rarity of solids 
in nature having spontaneous currents in their ground 
staters. 



A. Summary 

The central results of this paper are as follows. We 
emphasize first that our microscopic Hamiltonian was 
limited to (mainly spinless) models with interactions 
and one- fermion hopping terms. We did not explore 
the possibilities of correlated hopping, which were al- 
ready known to be conducive to the "d-density-wave" 
current orderi^*^. Ultimately such terms come micro- 
scopically from higher-order processes in fermion hops; 
thus, within our picture, related terms might be accessed 
by expanding our canonical transformation (Sec. lHll and 
Appendix!^ to higher orders, producing effective inter- 
actions with higher powers of pseudospin. 

We showed that by tuning the parameters and the 
underlying geometry of a toy spinless Hamiltonian, we 
can make a system acquire spontaneous currents, charge 
or bond order. The crucial ingredient of our analysis 
was the pseudospin mapping (Sec III Bp and degenerate 
second-order perturbation theory ( Appendix [X)) . We saw 
that bond ordering is naturally disfavored in ladders (Sec 
IHI Bp : while, for tube and square lattice, the ground 
state can acquire possibly coexisting charge, bond or cur- 
rent order in a symmetry breaking fashion (Sec IIII Al and 
IHI Bp . Since, the tube is quasi one dimensional, the cor- 
relation will actually be power laws; but for the square 
lattice, we will have true long range order. Furthermore, 
the pattern of currents corresponding to the spontaneous 
currents carrying ground state is same as the d-density 
wave state (Sec. IV Ap . The Nersesyan map provided a 
way to extend our strong-coupling result to that of an 
uniform case for the ladder (Sec IV B|) . 



B. Relation to three-band model and other real 
systems 

The current contact of our topic with real systems is 
in the three-band model of cuprates^i^^. A recent paper- 
argued (by mapping to a 2-channel Luttinger liquid and 
then analytic perturbation) that a ladder version of the 
three-band model of cuprates has long-range order with 
a current pattern similar to Varma's state. This claim 
was brought into question by a subsequent DMRG cal- 
culation^ on the same model: the current-current cor- 
relations were seen (numerically) to decay with a power 
law. But that, of course, indicates the presence of gap- 
less excitations, like the Goldstone mode of a continuous 
symmetry; it would not expected for the Ising-like sym- 
metry of a (pseudo)current operator, unless the system 
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is accidentally at a critical point. 

In our spinless model, it was easy enough to get current 
order degenerate with bond-density order, or to stabilize 
the latter, but quite hard to stabilize just current order. 
Could this be going in the ladder system of Refs. and 
|4^ . in which the current operator is just one component of 
a larger order parameter with XY symmetry? We warn, 
however, that our results also suggest that ladders (ap- 
plied as an approach to square lattice models^l) are plain 
deceptive. The key terms stabilizing currents depended 
on the very anisotropy of the ladder. 

For that reason, it is interesting that Ref. 8 do 
claim to obtain currents in a variational and truly two- 
dimensional calculation. So far, there is no clear picture 
of which interactions are crucial to causing the order. A 
strong-coupling caricature of the three-band lattice (or 
ladder) in the spirit of our models might illuminate this 
point. 
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method and the results that are used. Say that H.^ rep- 
resents the strong-coupling limit, which is assumed to 
have eigenenergies {Ea} with a large ground state de- 
generacy, which however is split by a small perturba- 
tion 7i' (that has no matrix elements between degener- 
ate states of 7i°). We desire to project our problem onto 
the ground-state ("zero") subspace of 7i". The usual 
way to accomplish this is canonical transformation: let 
H = e^^He"^^ , where we require Ha/3 = for any ma- 
trix element connecting the zero subspace to other states; 
then, we can restrict our Hilbert space to the span of 
e^'^ja), where \a) was one of the ground states. 

To lowest order in Ti', the standard canonical trans- 
formation is given by Sap = when states a and P are 
degenerate, and Sap = iTL'^p/ {Ep — Ea) otherwise. Then 
the effective Hamiltonian is given by = H + SH, where 

(Al) 

is the off- diagonal second-order perturbation correction. 
Similarly, an operator O is transformed to O = O + SO, 
where (to lowest order) 



APPENDIX A: CANONICAL 
TRANSFORMATION 

Since we use the method of Canonical Transformation 
for our calculations, here is a very brief summary of the 



SOap = i 



xp 



OaxH'y 



A7 



(A2) 



^\Ex^Ea Ex-Ep, 
For our problem, 7i''°P + is the strong coupling 
limit and TL' is the small perturbation. 
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